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0^ ■ Abstract. In this contribution we discuss gravitational effects of global 

scalar fields and, especially, of global topological defects. 
We first give an introduction to the dynamics of global fields and the 
P ! formation of defects. Next we investigate the induced gravitational 

'^ I fields, first in a fiat background and then in the expanding universe. In 

fiat space, we explicitly calculate the gravitational fields of exact global 
monopole and global texture solutions and discuss the motion of photons 

^ ■ and massive particles in these geometries. We also show that slowly mov- 

ing particles and the energy of photons are not affected in static scalar 
field configurations with vanishing potential energy. In expanding space, 
we explore the possibility that global topological defects from a phase 
transition in the very early universe may have seeded inhomogeneities in 
the energy distribution which yielded the observed large scale structure 
in the Universe, the sheets of galaxies, clusters, voids ... . We outline 
numerical simulations which have been performed to tackle this problem 
and briefly discuss their results. 
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1. INTRODUCTION 

During this meeting we have learned about phase transitions, the formation of topo- 
logical defects during phase transitions, and the dynamics which usually leads to 
certain scaling laws for the density of defects and their correlation length. In most 
of the previous talks gravitation has been disregarded; it was unimportant for the 
examples under consideration. In these lectures we want to discuss the gravitational 
interaction of global defects with matter and radiation. 

As we shall see, the gravitational coupling strength of topological defects is of 
the order of 

^ = GT^ = {Tjm,,f , (1) 

where T^ is the symmetry breaking temperature and nipi = 1/vG ~ lO^^GeV ~ 
10^^ K is the Planck mass {Ti = c = kBoitzmann = 1 throughout). For gravitation to 
become important, the symmetry breaking scale thus cannot be too far below the 
Planck scale. In the electroweak phase transition, e.g., with T^ ~ lOOGeV, fi ~ 10^^"^ 
and gravity can be ignored completely. 

We shall see later, that topological defects might be responsible for cosmological 
structure formation, if they form during a phase transition at Tc ~ lO^^GeV , n ~ 
10^®. This energy coincides roughly with GUT scale. Certainly, this energy scale can 
never be probed directly by accelerators or any present day astrophysical events like 
supernovae. There are thus justified doubts if we will ever have a detailed picture of 
the physics taking place at these energies. On the other hand, if the ideas explored 
here turn out to be correct, and topological defects due to a phase transition in 
the very early universe have triggered cosmological structure formation, then the 
galaxy distribution in the universe and the anisotropies in the cosmic microwave 
background (CMB) may be relics of physics at GUT scales! 

Since the generation and evolution of defects is quite a generic feature, we 
hope that the main results will not depend very sensitively on the detailed physical 
model. In that sense, I think, the scenarios discussed later in these lectures should 
be regarded as a kind of toy models which we hope are capable of capturing the main 
features, but we should not expect them to make predictions / verify observations 
to much better than within a factor of two. 

In the next section we introduce some generalities on the dynamics of global 
fields and defect formation. We mention some important results from homotopy 
theory, present the a-model approximation for the dynamics and discuss Derrick's 
theorem. In Section 3, we study gravitational effects of global fields in fiat spacetime. 
We calculate the gravitational influence of a global monopole and a global texture 
on matter and radiation. In Section 4, we investigate global defects in expanding 
space and especially the possibility that they may seed the formation of large scale 
structure in the universe. We shortly discuss cosmological perturbation theory, the 
Harrison Zel'dovich spectrum and numerical simulations of structure formation. We 
conclude in Section 5. 



2. GLOBAL FIELD DYNAMICS AND DEFECT FORMATION 

We consider a scalar field (order parameter) </> with Lagrangian density 

c = ]^d,<p-d^<p-v{4>), (2) 

G V, where V is a finite dimensional vector space and ■ denotes a scalar product 
in V. Here "scalar field" does not refer to the number of components of the field 
but to the transformation of under rotations of physical space: Be i? a rotation, 
then (L)(i?)0)(x,t) = (j){R-^x,t). 

If we quantize the field at finite temperature, we can take into account the 
interactions of the particles with the thermal bath by replacing V by an effective 
potential Vr- The precise form of Vr depends on V and on the interactions of 
with other particles, fermions, gauge bosons ... . In the simplest situation with only 
the scalar field G R^ and 

V = ^A(0^ - vr , (3) 

one finds in one loop approximation the high temperature corrections 

Vt = ^A(0^ - r^r - n^T^ + M^^!^T^ + 0(T) + V.r^o) , (4) 

where n denotes the number of helicity states of the field. At high temperatures, 
T^ '^iffi the only minimum of Vt is the field value = 0. As the temperature drops 
below the critical temperature T^ = 27], additional minima at (0)^ = rf{l — T"^ /T"^) 
develop and the vacuum manifold Mt (space of minima of Vt) becomes an (A^ — 1)^ 
sphere. 

More generally, we assume C to be invariant under the action of some compact 
Lie group G on V, which leaves only = invariant. If the vacuum manifold 
consists only of the invariant element = 0, the symmetry is unbroken. Since Vt 
is temperature dependent, at some other temperature, the vacuum manifold, A'r 
may become non trivial and contain an element 0o 7^ 0. Since Vt is invariant under 
G, the whole orbit {g(t)o\g ^ G} then belongs to A'r. If the symmetry group G is 
maximal, Mt consists just of the orbit of 0o which is given by G/ H, where H G G 
denotes the invariance group of 0o, H = {h G G\h(j)o = 0o}- The symmetry G is 
then spontaneously broken to the remaining symmetry group H. 

Even though the opposite case can also occur [0], we shall generally assume 
that the symmetry is restored at high temperatures, T > Tc and spontaneously 
broken at lower temperatures T < T^. If the temperature then falls below the 
critical temperature, T^, and if N't is topologically non trivial, defects of dimension 
d in spacetime can form via the Kibble mechanism [Q]. The collection of the different 
topological defects possible in four spacetime dimensions is presented in table |l[ 

Since the reader is probably quite familiar with the appearance of domain 
walls, strings and monopoles, let me just briefly explain textures: We consider a 



Homotopy 7r„, dimension in spacetinie= d=4-l-n 


appearance 


vro(Ar) ^ 

M disconnected 


walls form 


d = ?> 


sheets in space 


TTi (A/") 7^ A/" contains 
non shrinkable circles 


strings form 


d = 2 


lines in space 


7r2(A/') 7^ A/" contains 
non shrinkable 2-spheres 


monopoles form 


d=l 


points in space 


TTslA/") 7^ A/" contains 
non shrinkable 3-spheres 


textures form 


d = 


events in spacetime 



Table 1. Topological defects in four dimensional spacetime 

field configuration, which is asymptotically constant (as it has to be if we require 
the field to have finite energy). At fixed time t, cj) can then be regarded as map from 
compactified space R^ = R^ U {oo}, with 



(oo) 



lim ( 

|x|— >oo 



(x). 



Since R^ is topologically equivalent to S^, we can now regard a vacuum configuration, 
{/), as a map from S^ into M . If the image, 0(S'^) (which is of course topologically 
again S^) is not contractable in A/", the configuration cannot evolve into the trivial 
one, =constant, without leaving the vacuum manifold. Such a configuration is 
called texture. If has finite energy. Derricks theorem tells us that it will shrink 
and eventually evolve into the trivial configuration, leaving the vacuum manifold at 
some spacetime event (with extension of the order the inverse symmetry breaking 
scale). This event is the texture singularity. 

We now want to state a few facts from homotopy theory which are commonly 
used throughout. Proofs and further information can be found in [^, ^, ^. 

In contrast to homology groups, there exists no general algorithm to calculate 
homotopy groups 7r„ for n > 1. Although a lot of research has been carried out, not 
even all homotopy groups of the sphere are known! 

The following results from homotopy theory of Lie groups are often useful: 
Since every compact connected Lie group is a product of some of the following 
groups: 

U{n), SO{n), SU{n), Sp{n), Spin{n),G2,F4,EQ,E7 or Es , 

it is sufficient to discuss these groups. Here Spin) denotes the symplectic group of 
dimension 2n, Spin{n) denotes the spin group of dimension n, i.e., the universal 
covering group of SO{n) and G2 ... Eg are the exception groups. Be G one of the 
above groups, then 

7ri{U{n)) = Z, 7ri(50(n)) = Z2 and 7ri(G') = for all others. 



7l2{G) = 



7r3(G) = Z, if 6-^50(4) and tx^{SO{A)) = Z ®Z . 
TikiUin)) = for all A; > 1 . 

Corresponding identities for direct products follow from 

vrfc(Gi X G2) = TTkiGi) © TTkiGi) . 

The main tool to determine 7r„ are exact sequences. 

Definition: Be A^ a sequence of sets and (pn a mapping from An to A„+i. 
The sequence 

... /in -^ ^n+1 -^ ^n+2 ••• 

is called exact if 

ker{ipn+i) = im{ipn). 

Theorem: For a subgroup H G G the sequence 

...7ln{H) -^ 7r„(G) A Tln{G/H) ^ 7r„_i(/J) ^ 7r„_i(G)... 

is exact. Here j and i denote the trivial inclusion map and d is the boundary map. 
Since vr2(G') = 0, we thus obtain 

7!'2{G/H) = 7ri(iJ) for all simply connected groups G; 

i.e. for all groups with 7ri(G') = 0. 

2.1. Tlie a~ model approximation 

If the system under consideration is at a temperature T much below the critical 
temperature, T -C T^, it becomes more and more improbable for the field to leave 
the vacuum manifold. will leave the vacuum manifold only if it would otherwise 
be forced to gradients of order (V0)^ ~ XcfP'rf, thus only over length scales of order 
/ = \/{-\/\ri) = m~,^ {I is the transversal extension of the defects). If we are willing 
to loose the information of the precise field configuration over these tiny regions (for 
GUT scale phase transitions / ~ 10~^°cm as compared to cosmic distances of the 
order of lMpc~ 10^'^cm !!) it seems well justified to fix to the vacuum manifold 
N . Instead of discussing the field equation from (|^), 

dV 
□*+g^ = 0. (5) 

we require G M. A/" C V is a Riemannian submanifold with the induced scalar 
product. The remaining field equation □(/) = then just demands that 

d) : M ^U 



is a harmonic map from spacetime Ai into Af. There exists a waste mathematical 
hterature on harmonic maps and their singularities which might be useful for us and 
should be explored 0. 

The topological defects we are interested in are singularities of these maps. 
When the gradients of become very large, like, e.g., towards the center of a global 
monopole, the field leaves the vacuum manifold and assumes non vanishing potential 
energy. If </> G A/" is enforced, a singularity develops by topological reasons. 

In the physics literature harmonic maps are known as a-models. They were 
originally introduced because of their similarities with non Abelian gauge theories 
(the corresponding field equations also contain non-linear gradient terms). The 
action of a a-model is given by 

S.= f g'"'d,(P^d,(P''hAB{(p)J\g\d'x , (6) 

where Hab denotes the metric on A/" and g^^ is the metric of spacetime. Let us 
consider once more the 0{N) example: 

V = l\i<p' - v'f , 0eR^. 

We now fix (p to lay in the vacuum manifold, S^^^ with radius r], by introducing a 
Lagrange multiplier. 

C, = d^<P ■ a^0 - a(02 - r^^) . 

Variation w.r.t </> yields 

□0 + 2a0 = . (7) 

We multiply (^ with to obtain a = —(f) ■ □0/(2?7^). Inserting this in ^ we obtain 
the field equation 

□^_(^^.0. (8) 

In other words, the projection of D^ onto the hyperplane tangent to the sphere has 
to vanish, i.e., D^ = on Af. In terms of the dimensionless variable f3 = (p/r], (|^) 
reads 

DP-{P.DP)P = 0, (9) 

which shows that the cr-model is scale free. 



2.1.1. Analytic flat space solutions 

A global string along the z— axis is described by the field configuration e R^: 



(f){x,y,z) = rjep = r]{cosip,smip), cr-model (10) 

(l){x,y,z) = fs{p)ri{cosip,siinp), full field eqn. (11) 

where fs satisfies 

fs + -fs'\fs + ^{f's-l)fs = 0, (12) 

p p z 

with boundary conditions /^(O) = and fs{oo) = 1. Here p is the cylindrical radius 
and ip the polar angle, {x,y) = p{cos(p,sm(p). 

A spherically symmetric, static global monopole is described by the field 
configuration (f) G R'^ with 



(f){x,y,z) = rjBr = ri{sm 9 cos ip, sin 9 simp, cos 9), a-model (13) 

(f){x,y,z) = f M {r )ri {sin 9 cos ip, sin 9 simp, cos 9), full field eqn. (14) 

where /m satisfies 

I'm + -I'm - -Jm + ^{11 - ^)fM = , (15) 

r r Z 

with boundary conditions /a/(0) = and /a/(oo) = 1. The equations for fs and /m 
can only be solved numerically. 

A spherically symmetric global texture is described by the field configuration 
e R"^ with 

= ?7(sin X sin 9 cos ip, sin x sin 9 sin ip, sin x cos 9, cos x) ■ (16) 

With the ansatz x = xij^i)) the cr-model field equation (^ leads to 

Since the a-model is scale invariant, we further require x = xiv) with y = t/r. In 
terms of this scaling variable, the equation of motion for x reduces to the ordinary 
differential equation 

{y' - l)x" = sin(2x) , (17) 

with exact solutions 

Xiy) = 2arctg(±y) ± nn . 



Figure 1. The function x is shown for t < 0, sohd curve, and for t > 0, dashed curve. 
X goes from to vr for negative times, i.e. the configuration winds once around S^, 
and from vr to |7r and back to tt for positive times, i.e., no winding, 
^paragraph. 

These solutions were originally found by Turok and Spergel |]^ . To obtain a solution 
which winds around the three sphere at negative times and collapses at t = 0, we 
patch together x as follows: 

vfv/l = / 2arctg(|/) + TT , -oo < y < 1 

^^^' \ 2arctg(l/i/) + TT , 1 < y < oo ^ ' 

The behavior of x as function of r for positive and negative times is shown in 
Fig. 1. The kink at r = t for positive times is due to the singularity of the a-model 
solution at r = t = 0. It would be softened in a solution of the full field equations. 
Physically, this kink represents the wake of Goldstone bosons in which the massive 
mode at r = t = has decayed and which now travels out with the speed of light. 
One easily sees that the energy of these three configurations diverges. For a large 
ball of radius R one finds 

Estring{R) = JT^,{z,p) = ]^j{V<t>f OC i?log(i?r/) , (19) 

Emonopole{R)= JT^{r) = ]^jiy(Pf (X R, (20) 

Ete.ture{R) = jTl^{r) = ^J[^^ + {V<j)f] (X R for i? » t . (21) 

Before we go on to discuss the gravitational effects of these solutions, let me 
briefly note some thoughts concerning Derrick's theorem. 



2.2. Derrick's theorem 

Since it is so simple and beautiful, let me state the theorem with proof. [H 
Theorem: In d = 3 dimensions there are no non trivial static finite energy 
solutions for a scalar field whose potential energy is bounded from below. 

Proof: For static configurations, the variation of the action can be replaced 
by the variation of the energy, E. 



with 



E = J[^{V4>f + V{4>)]£x = h + h , 
h = j]^{V<pfd?x , l2 = Jv{<f))d'x 



Without loss of generality, we may assume V^ > (otherwise, consider E — Knin)- 
Then Ji > and I2 > 0. We assume now (f){x) be a non trivial solution and consider 

For the scaled configuration we have 

/2(0a) = X^h . 



the scaled configuration (f)x{x) = 


(j){X-^x 


h{x) -= 


= /i(0a) = A/i and 


hix) = 


Therefore 






dxEx- 


=1 = Ji + 3J2 > . 





This contradicts our assumption of beeing a solution. □ 

From this we can immediately conclude that our solutions for global strings and 

monopoles discussed before must have infinite energy. But also the time dependent 



texture solution has infinite energy (pl|). 

Perivolaropoulos0 has put forward the following argument: In the cosmo- 
logical context we should truncate the energy at some large radius R, the horizon 
distance or the distance to the next defect. Then the variation of the scaled energy 
yields 

dxE\x=i =h + 3/2 - Rdnih + h) , 

which, due to the negative term, can vanish. The second variation shows that a con- 
figuration with vanishing first variation does represent a minimum of the truncated 
energy and thus is stable against shrinking and expansion. 

But of course this argument does not explain the existence of the string and 
monopole solutions considered previously. Furthermore, the argument would also 
allow for stable static texture solution (with infinite energy). There have been 
some analytical and numerical arguments [|TU], ^, 0, |T^, [I^ , that it is the winding 



condition that renders the textures unstable. For winding number n > 0.5 textures 
tend to shrink and for n < 0.5 they tend to expand. Nevertheless, in my opinion, a 
clear understanding of the numerical finding that there exist stable static (infinite 
energy) string and monopole solutions, but probably no stable static texture solution 
is still missing. 



-^ flU ^MV ^UY 


^UlA.u^AH' ^ V • 


We set 






= \{^' + iy<pf) 


P=il? 


= ^(0^-wn 


^ij -^ ij 9ijP 


= di(j)dj(f) - -g,j{S/4>f . 



3. GRAVITATIONAL EFFECTS OF SCALAR FIELDS IN FLAT 
SPACETIME 

3.1. Generalities 

The energy momentum tensor of a scalar field configuration in the a-model approx- 
imation is given by 

T(f = d,ct> ■ 5.0 - ]:g,udx<P ■ 9> . (22) 



(23) 
(24) 

(25) 

For static global field configurations p + 3p = 0. This indicates that static 
global field configurations, like an infinite straight string or a hedgehog monopole, 
do not gravitationally attract nonrelativistic particles. 

To discuss the gravitational effects of test particles and radiation in general, 
we have in principle to solve Einsteins equations, 

G,, = MGTl:^ , (26) 

and investigate the geodesies in the resulting geometry. For a typical field coherence 
length Z, we have SnGT^,^ ~ SttG?]^//^. For a GUT phase transition this is of the 
order of 10~^//^ — 10"'^//^. The induced changes of the metric will thus be small, of 
order 10~^ — 10~^, and we can treat gravity in first order perturbation theory. I.e., 
we insert in eqn. ( ^B]) the unperturbed, fiat spacetime, energy momentum tensor 
and equate it to the Einstein tensor G^u obtained from first order corrections to 
the flat metric (or, in the cosmological context to the Friedmann Robertson Walker 
metric). 

3.2. Spherically symmetric field configurations 

For the sake of simplicity, we now restrict ourselves to spherically symmetric con- 
figurations. In first order perturbation theory, the metric can then be parametrized 

by 

g = -{l- 2^)dt^ + (1 - 2<!?)6ijdx'dx^ . (27) 

The linearized Einstein equations yield 

-A<l> = AnGp (28) 

-A(<l>-$) = SvrGAn , where didjU - -AU = iiij . (29) 



(In the spherically symmetric case it is always possible to find such an anisotropy 
potential 11.) For ordinary matter, p ^ ttj^, and thus $ = — \E'. ^ is the relativistic 
analog to the Newtonian gravitational potential, and slowly moving matter only 
couples to ^ . Using the equation of motion (§) for 0, one can show that for static 
configurations 

An=i(V0)^ = ip. (30) 

Eqn. (^) and (^) then yield \E' = 0. This shows again that nonrelativistic matter 
is not affected by static global field configurations. 

It is easy to calculate the connection coefficients (Christoffel symbols) from 
Ansatz (^Tf ). Inserting them into the geodesic equation for a photon moving with 
four velocity 

n = {l,n) + Sn , p = En , 

one obtains in first order perturbation theory 

6n^ = -2-^\{ + f {^ - ^)dX (31) 



Sni = I di{^ - ^)rfA , (32) 

where the integrals are performed along the unperturbed photon trajectory. The 
meaning of these quantities is the following: We consider an emitter/ observer of a 
light ray moving according to the velocity field 

M = (1 + ^, v) with v^ < 1 . 

The 0-component of u is determined by the condition v^ = —\. The energy shift of a 
photon relative to emitter and observer is generally given by 5E = {p-u){f) — {p-u){i). 
In our situation this yields 

6E/E = 6{u-n) = n- v|f + *|f - / ($ - 4/)dX . (33) 

Ji 



The first term on the right hand side of (|33D is the usual, special relativistic Doppler 



term. The second term is due to the difference of the gravitational potential at 
the position of the emitter and observer, and the third term is a path dependent 
contribution due to the change of the gravitational potentials during the passage 
of the photons. Since for static scalar fields \I' = \I' = $ = 0, the gravitational 
contributions to 6E vanish in the static situation. 

Eqn. (^) is related to light defiection. Be e the radial unit vector. The 
defiection of a light ray emitted at position i and observed at / is then given by 



a 



6{n ■ e){t) - 6{n ■ e)(/) = - [^ 9,(<l> - ^)e,dX . (34) 

Ji 



For the gravitational field from ordinary matter ($ = — \E'), we recover the old result 
by Einstein (with the correct factor of 2). 

For a slowly moving massive particle in a weak spherically symmetric gravita- 
tional field, we make the ansatz 

u = (1,0) + Su . (35) 

From the geodesic equation, we then obtain 

5u' = - j dS!dt . (36) 



Since \E' vanishes in the static case, we find that slowly moving particles are not 
affected by static field configurations. 

Taking into account also (|55|), we thus have proven the following 
Theorem: In static scalar field configurations with negligible potential energy, the 
gravitational redshift of photons and the gravitational acceleration of slowly moving 
particles vanish. 

The gravitational field of static configurations thus affects matter and radiation 
only by deflection of relativistic particles. 

3.3. Two examples 

To be somewhat more specific, we now want to insert into ( p3|) and ([3^ ) the global 
monopole and global texture solutions obtained in the last section. 

Global monopole: From the linearized Einstein equations we find for the 



hedgehog monopole solution |T5| 



^ = , $ = -SvrGr/^ ln(r/Z) = -eln(r//) , (37) 

where we have set e = SixGr]'^ and / is an arbitrary constant of integration. We 
consider a light ray passing the monopole with impact parameter b. Its unperturbed 
trajectory is given by x(A) = An + 6e. Since the configuration is static, 6E vanishes. 
Inserting (|37D in ( |3^ ) yields the deflection angle 



aM = - I d^^e'dX = t I p^5?^^ = ^^ • (38) 



This result was originally found by different methods by Barriola and Vilenkin |T^ 
Global Texture: For the texture solution (p!8|) we obtain [l^ 



*4-m' --i-m- (-, 



For a light ray passing the texture with impact parameter b at impact time r {t 
r + A, r^ = 6^ + A^), we find the defiection angle 



ff 
ar = - / ($ - ^),ie*(iA 



r 

J — ( 



^ dX 



62 + 2A2 + 2Ar + t^ 
b 



en 



VWT^ 



(40) 



This result was first obtained in |T^ 



To calculate the energy shift of a photon passing the texture, we have to 'renor- 
malize' the result obtained from naively inserting (|39D in (^). Due to the unphysical 
infinite energy of solution (|18D , the energy shift contains a divergent logarithmic term 



which we neglect in the ~ sign in eqn. (|41|). (In [|T^ this renormalization is discussed 
in some detail.) 

^.>Hf + /(*_*),A«..^=L=. (41) 

This result was first obtained in |^. 

The interesting difference between the results for monopoles and texture is 
due to the time dependence of the latter. This, first of all, yields a non vanishing 
energy shift for the texture. An observer receiving photons from behind a collapsing 
texture sees them first redshifted (if they pass the texture before collapse) and then 
blueshifted (see Fig. 2). An observer in perfect alignment with a background quasar 
and a global monopole sees the quasar image as Einstein ring with fixed opening 
angle. In the case of a global texture, the Einstein ring opens up some time before 
texture collapse, reaches a maximum opening angle of the same order of magnitude 
as in the monopole case and then shrinks back to a point [1T7| , |15|. 

The gravitational field of our global texture solution (|T8|) also accelerates slowly 
moving particles. Inserting (|39| ) in ( ^6]) leads to the wellknown result 0, ^ 



v(/) - v(z) = -evrer . (42) 

Slowly moving particles around a collapsing texture thus acquire a net infall velocity 
of amplitude en. 



4. GLOBAL DEFECTS AS SEEDS FOR COSMOLOGICAL STRUC- 
TURE FORMATION 

Many observational results, like Hubble expansion, primordial nucleosynthesis, the 
isotropy and the thermal spectrum of the cosmic microwave background, confirm 
the idea that on large scales the Universe is homogeneous and isotropic. On large 



Figure 2. The temperature fluctuation, AT/T, induced by a spherically symmetric 
collapsing texture as function of the impact time of the observer. The solid line 
shows the result in expanding space, the dashed line is the fiat space result. The 
collapse time of the texture is tc ~ 20 (in arbitrary units). The difference of the 
two curves is due to the existence of horizons in expanding space: Photons, which 
pass the texture long before of after the collapse are not influenced in expanding 
spacetime, but acquire the maximum energy shift in flat spacetime. 
^paragraph. 

scales, the observable Universe is thus well approximated by a Friedmann universe, 
which evolved from a very hot thermal state, the big bang, by adiabatic expansion. 

On smaller scales, clearly, the Universe is lumpy. Laborious mapping of the 
3d galaxy distribution has shown that this dumpiness persists on scales up to (30 - 
50)/i"^Mpc. The galaxies themselves are arranged in relatively thin sheets surround- 
ing seemingly empty voids of diameters up to 50/i^^Mpc. {IMpc ~ 3.2 x IQ^ly ^ 
3.1 X lO^^cm) (see Fig. 3). 

With the help of the Cosmic Background Explorer (COBE) satellite, anisotro- 
pies have been found also in the cosmic microwave background (CMB) which are on 



the level fH 



VAT\n _5 



A ( f -^ ) ){0) ^ 10"^ , on all angular scales 6* > 7° 



These findings support the old idea of Lifshitz [21] that the cosmic structures might 
have formed by gravitational instability from small initial fluctuations. 

Cosmological perturbation theory shows, that perturbations in the radiation 
field can not grow substantially. Therefore, AT/T yields the amplitude of initial 



Figure 3. The distribution of more than 1000 galaxies of the CfA catalog. A slice of 
the universe, about lO'OOOkm/s deep in redshift space and about 5° thick is shown 
(from Geller and Huchra ||19| ). 
^paragraph. 

fluctuations {6p/p)in ~ 3AT/T. On the other hand, perturbations in pressureless 
matter {p -C p, cosmic dust) grow roughly by a factor ao/cteg = -^eg + 1, where 
a denotes the scale factor of the universe, a subscript o denotes present time and 
eg denotes the time of equal matter and radiation density. If the matter content 
of the universe is given by baryons only, Zeg < 10^, and our naive estimates lead 
to perturbations which are roughly by a factor 10 too small to yield the observed 
structures. However, if we assume that the universe is dominated by dark matter 
leading to critical density, fi = 1, we have z^g ~ lO'', of the correct order of magnitude 
to lead to the nonlinear clustering observed today. 

There remains one basic ingredient to the gravitational instability picture: 
How did the small initial perturbations of order 10"^ — 10~^ emerge? Presently two 
mechanisms are primarily investigated: 

• Quantum fluctuations 'frozen in' as classical perturbations of the energy den- 
sity after an epoch of inflation. 

• Topological defects from a phase transition in the early universe. 

In this workshop, we concentrate on the second possibility. We have seen in the 
last section that topological defects yield gravitational perturbations of the order 
of SttGt]^ = e. To obtain e = 10~^-10~^, we need a GUT scale phase transition, 
r] ~ lO^^Gel/. 



4.1. Scaling 

Let us now assume that on large scales the Universe can be described by a Friedmann 
universe with vanishing spatial curvature, Q = 1. The metric of spacetime can then 
be given by 

ds^ = a^{-dt^ + 6ijdx'dx^) . (43) 

Here a is the cosmic scale factor and t is conformal time. It is related to the cosmic 
time, tcos) which has elapsed since the big bang by 

tcos{t) = / a{t')dt' 
Jo 

(see also contribution by T.W.B. Kibble). To be relevant for structure formation, 
topological defects must make up an approximately constant fraction of order e of 
the total energy density of the universe. In the cosmological context, we then say 
that the defects obey scaling. Let us estimate the energy density of global defects, 
neglecting the potential energy: 



'{p^ef) = h'^<p?) + hd,cf^r)-vVt' 



where we have assumed that changes typically over a horizon scale. On the other 
hand, from the Friedmann equation (assuming the scale factor a to obey a power 
law), we have 

a p = ^^{a/a) 



rsj 



so that {pdef)/p ~ e- 

From this result one might conclude that all global defects obey scaling. But 
the above argument is somewhat too simplistic as the case of global strings shows: 
Let us approximate the energy of a global string inside one horizon volume by the 
corresponding energy of a straight cosmic string in fiat space: 



E{t) = 271 f '"' dzrdr{rf/r^) = 2nri\os\og{t,osr]) 

Jq 



and thus a^{pstr) ~ rf / 1^ \og{atri) . 

In the case of strings, we thus obtain a logarithmic correction term which, for a GUT 
scale phase transition, amounts to a factor of approximately 150 today. For higher 
0{N) defects like monopoles, textures and 0{N) models with A^ > 4, the scaling 
behavior becomes clean (see Fig. 4). 

In the case of local defects, only cosmic strings obey scaling. Monopoles stop 
interacting soon after formation and then scale like massive particles: 



um ~ l/tl{ac/af , pm = mMUM oc a 



-3 



Figure 4. The scaling behavior of (p + 3p)a^ found numerically in (128)^ simulations 
is shown for four different 0{N) models. Time is given in units of the grid spacing 
Ax. For comparison, the dashed line oc 1/t^ is shown. For A^ > 3 scaling is very 
clean until t ~ 80, where finite size effects can become important, 
/.paragraph. 



The universe at GUT scale is radiation dominated, p (x. a ^. Therefore, soon after 
the phase transition pm ^ p leading to 



^M{to)h' ~ lO^\TjlO^^GeVy{mM/lO^^GeV) > n^h' ! 



This is the famous monopole problem in cosmology |^ . The reason why this repre- 
sents a serious problem is the following: Imagine some simple, compact grand unified 
group G, like SU{5), breaking (in one or several steps) to H = SU{3) x SU{2) x U{1). 
The existence of monopoles at the end is then determined by the exact sequence 

T^2{G) ^ ^2{G/H) ^ ^,{H) ^ 7^(6-) . (44) 

Since it2{G) = 0, and for a simple group also vri(G) = we find 

n2iG/H) = n,{H) = Z. 
Monopoles thus always form. By the analogous sequence for ni^G/H), 

= 7ri(G) ^ 7ri(G/iJ) ^ no{H) = , (45) 

we conclude that no strings form. The monopoles are thus not connected by strings 
and are stable. 

This is a beautiful example showing that observations of the present universe 
can lead to predictions about high energy physics and cosmology at GUT scale. The 
most simple GUT scenario is not compatible with standard cosmology. One either 
has to invoke a period of inflation or change the GUT idea [p^] . 



4.2. Cosmological perturbation theory 

So far, we have only seen that the orders of magnitude come out reasonable for 
structure formation with topological defects from a GUT scale phase transition. 
We would hke to obtain more precise results. We want to simulate the evolution 
of defects, calculate the gravitational fields they produce, which in turn affect the 
distribution of matter and radiation. We want to calculate the induced anisotropies 
in the cosmic radiation field and in the matter distribution, (AT/T)(to,x, n) and 
{5p/p){to,x), Vpec{to,x)- 

An important tool for this calculation is cosmological perturbation theory. We 
do not develop it here, but just mention the basic equations which determine our 
problem. For more details see, e.g., P3| . 

• The equation of motion for the scalar field: 

dV 
D(j) -\- —— = (potential model), or (46) 

0(p 

□0 - (p^—^-^ = (sigma model) . (47) 



• The perturbation of the energy momentum tensor: 

5T^. = T^. (</)) + <5T™«"^^. (48) 

• The hnearized Einstein equations: 

5T^, = 5G,, . (49) 

• The equations of motion hnearized about the Friedmann background: 

- The Liouville equation for photons 

p'^5^/ + r>V^ = o. (50) 

- The cold dark matter equation of motion, p = 0, 

T'^";.= 0. (51) 

For a perturbed Planck distribution, the Liouville equation can be cast into a per- 
turbation equation for the temperature only |jl5|, |2^: Be x the observer position and 
n the direction of observation. If we set T(x, n) = T(l +m(x, n)), the perturbation 
equation corresponding to (|50D can be expressed as 



{dt + n'di)x = -Sn'd'Ei, - n'ri^eindiBij , (52) 

where 

X = Am + (monopole term + dipole term) , 

and Eij, Bij denote the electric and magnetic parts of the Weyl tensor. If we are 
only interested in the spherical harmonic amplitudes aim of rn{n) for harmonics 
higher than the dipole, / > 2, it is thus sufficient to determine A~^x- Fo^' ^ fixed 
observer position a monopole term can not be distinguised from the background 
temperature and a dipole term can be attributed to the peculiar velocity of the 
observer. Therefore, monopole and dipole terms anyway do not contain information 
on the temperature fiuctuations. 

From ( |51| ) we obtain a perturbation equation for the energy density perturba- 
tions, D of the dark matter. 

D + {d/a)D - iTrGa'^pdniD = 47rG(p'^ . (53) 

From X, we can obtain m = 6T/T by inverse Laplacian. The first information 
to be compared with observations are the power spectra or, correspondingly, the 
auto-correlation functions of 6T/T and D. We expand 6T/T in spherical harmonics: 

((5T/T)(to,x, n) = ^ aim(x)yi„(n) . 

Lm 



The power spectrum of 5T /T is then given by 

where n^ is the number of observer positions x averaged over, and 2/+1 is the number 
of values —l<m<l. One easily finds the temperature correlation function p| 



(((5T/T)(n)(5T/r)(n'))„.„'=cose = ^Y.i'^l + l)ciPi{cose) . (55) 

Pi denotes the /th Legendre polynomial and () indicates averaging over positions 
and over all directions n, n' with relative angle 9. 

The power spectrum of dark matter perturbations (called 'structure function' 
in condensed matter physics) is the Fourier transform of the correlation function. 
Indicating Fourier transforms by a tilde, we have 

P{k) = \D{k)\'^ = C{k) , where (56) 

C(r) = (Z}(x)D(x + nr))x,„ 



is the correlation function [25 



4.3. The Harrison ZeI'dovich spectrum 

Let us assume that the only scale in the structure formation problem is the horizon 
scale. Then we expect the variance of the mass perturbation on this scale to be a 
constant. A, independent of time |2^, 



A = {\6M/M\')(^2./k=t) ~ k'\P{k,t = 27i/k)\ . (57) 

Once the perturbations 'enter the horizon', k > 27r/t, their behavior depends on the 
expansion law of the background spacetime. From a simple analysis of linear per- 
turbation theory one finds that perturbations cannot grow if spacetime is radiation 
dominated (Mezaros effect, [^), and they grow proportional to the scale factor a 
if spacetime is matter dominated. Let us denote by teq, aeq the conformal time and 
scale factor of the universe at the time when the energy density of radiation equals 
that of matter. During the matter dominated regime the scale factor grows like 
a{t) ex t^. Defining at = a{t = 27T/k), we obtain on scales which are subhorizon 
today (/c > 27r/to) 

~ / Ak-^/\ao/ak) = Ak'/\to/2ny , k < 2n/t,q 

^[l^,to) ~ I Ak-^/\ao/a,q) = Ak-y^z^q , k > 2n/t,q . ^^^^ 

The the Harrison ZeI'dovich spectrum can thus be approximated roughly by the 
form 



Figure 5. The points are the IRAS redshift space spectrum with f2 = 1. The box 
indicates the power spectrum inferred from the COBE DMR resuhs with spectral 
index n = 1. The sohd hue is the spectrum of a standard CDM scenario with Q = 1, 
normahzed to the real space variance of IRAS galaxies, as = 0.7 (this figure is taken 



from Fisher et al. [p9|]). 
^paragraph. 

with keg = 271 /teg. From large scale structure observations, the A;~^ behavior of the 
spectrum on small scales is approximately confirmed (the deviations on the smallest 
scales are probably due to nonlinear clustering). The bending of the spectrum on 



large scales is not yet observationally confirmed, see Fig. 5, [p^ . 



Correspondingly, one can show that for a scale invariant spectrum (0), the 



microwave background fluctuations behave like p5 



Numerical simulations and analytical arguments show that structure formation by 
global topological defects leads to an approximately scale invariant spectrum of 
perturbations. 

4.4. Numerical Simulations 



4.4.1. The scalar field: The quation of motion of a scalar field in expanding space 
is given by 

9'"'V,V.(l) + ^ = 0. (61) 

Defining P = (p/r] and m = v\r], ( |6TD yields for our 0{N) models in a Friedmann 
universe 

d^p + 2{d/a)dtp - V^j3 = ^a^m\p^ - l)/3 . (62) 

This equation as it stands is not tractable numerically in the regime which is in- 
teresting for large scale structure formation. The two scales in the problem are the 
horizon scale t and the inverse symmetry breaking scale, the comoving scale {am)~^. 
At recombination, e.g., these scales differ by a factor of about 10^'^ and can thus not 
both be resolved numerically. 

There are two approximations to treat the scalar field numerically. As we shall 
see, they are complementary and thus the fact that both approximations agree with 
each other within about 10% is reassuring. The first possibility is to replace {am)~^ 
by w, the smallest scale which can be resolved in a given simulation, typically twice 
the grid spacing, w ~ 2Ax. The time dependence of {am)~^ which results in a 
steepening of the potential is mimiced by an additional damping term, 2 (a/a) — *• 
ja/a, with 7 ~ 3 |^|. Numerical tests have shown, that this procedure, which 



usually is implemented by a modified staggered leap frog scheme |]3T|, is not very 
sensitive on the values of 7 and w chosen. With this method we have replaced the 
growing comoving mass am by the largest mass which our code can resolve. For a, 
say (256)^ grid which simulates the evolution of the scalar field until today, we have 
256 Aa; ~ to ~ 4 x 10-^ ''sec, so that w ~ 2 x lO-'^^sec, i.e., am ~ rj/zrec ~ lO^^GeV is 
replaced by about w~^ = 10~'^'^GeV\ 

We believe this mimics the behavior of the field, since the actual mass of the 
scalar field is irrelevant as long as it is much larger than the typical kinetic and 
gradient energies associated with the field which are of the order the inverse horizon 
scale. Therefore, as soon as the horizon scale is substantially larger than Ax, the 
code should mimic the true field evolution on scales larger then w. But, to my 
knowledge, there exists no rigorous mathematical approximation scheme leading to 
the above treatment of the scalar field which would then also yield the optimal choice 
for 7. 

Alternatively, we can treat the scalar field in the a-model approximation. This 
approach is opposite to the one outlined above in which the scalar field mass is much 
too small, since the a-model corresponds to setting the scalar field mass infinity. 

The cr-model equation of motion cannot be treated numerically with a leap 
frog scheme, since it contains non-linear time derivatives. In this case, a second 
order accurate integration scheme has been developed by varying the discretized 



action with respect to the field |^ . 



Initially, the field itself and/or the velocities are laid down randomly on 
the grid points. The initial time, t^ is chosen to be the grid size, t^ = Ax, so that 
the field at different grid points should not be correlated. The configuration is then 
evolved in time with one of the approximation schemes discussed above. 

The two different approaches have been extensively tested, and good agreement 
has been found on scales larger than about 2-3 grid sizes |Q Q. This is very 
encouraging, especially since the two treatments are complementary: In the a- 
model, we let the scalar field mass m go to infinity. In the potential approach, we 
replace it by ~ l/Ax ~ l/Uec ~ 1/100/y ~ 10~34QgY_ 

The integration of the scalar field equation is numerically the hardest part of 
the problem, since it involves the solution of a nonlinear partial differential equation. 

4.4.2. The gravitational perturbations: Once /?(x, t) is known, we can calcu- 
late the energy momentum tensor 

Tl£> = v'[{d,f^ ■ 9,(5) - \g,u{d^P ■ d^^ ■ (63) 



From eqn. (|53D we can further determine the perturbation of the dark matter en- 
ergy momentum tensor. The perturbed Einstein equations then yield an algebraic 
equation for the electric part of the Weyl tensor, Eij, and an equation of motion for 
the magnetic part of the Weyl tensor, Bij (see ||2^). The magnetic contributions. 



which consist of vector and tensor perturbations only, usually amount to about (10 
- 20) % of the electric contributions which are a combination of scalar and vector 
perturbations of the gravitational field, see Fig. 6. 

4.4.3. The perturbations of the cosmic background radiation and the dark 
matter: Using eqn. (|5^), we finally obtain 

ft 

From eqn. (|53|) , we can calculate D and the power spectrum P{k) = |-D(fc)p. 
Patching together simulations with different physical grid size, we can enlarge the 
range of comoving wave numbers k covered, see Fig. 7. 
By the decomposition into spherical harmonics, 

a^„^(x) = / i5T/T)ito,x,n)Yi^iii)dn , (64) 



f\n'A-\d^E,,) + n'n^A-\e,kid'B,,)}{t', x - n(t - t'), n)dt' 



and (plD, we determine the q's. Since the spectrum is close to scale invariant, it is 
entirely determined by the quadrupole moment, Q, 

An 
Q = —C2^a€ = QcoBE = (0.6 ± 0.1)10-^ . (65) 



Figure 6. The amplitude of the electric and magnetic source terms to the photon 
equation if motion are shown as a function of wavenumber k in arbitrary scale. For 
small wavelength (large scales) the magnetic part contributes about 1/4, decaying 
to roughly 1/10 on small scales. The quantities graphed are E = ^ Y.iid-' EijY , B = 
5 Hiji^jikd^BkiY . B is represented by the sohd hue. 
/,paragraph. 



Figure 7. The final dark matter spectrum of density fluctuations for 3 texture 

simulations with different physical grid sizes patched together. 

/.paragraph. 



The value of a above can be obtained by numerical simulations |35, ^, 32, 23| and 
is typically of the order 0.1 < a < 1. For textures one finds [^ 

e = %TxGif = (2.2 ± 1) X 10-^ . (66) 

The shape of the dark matter spectrum is again approximately determined by scale 
invariance. The integral of the dark matter perturbation spectrum over scales larger 
than R, determines the mass variation, cr^(i?) over these scales. The comparison of 
this dark matter mass variation with the observed variation of the galaxy distribution 
yields a scale dependent bias factor, b{R). The bias factors obtained this way are of 
the order 6 ~ 2 - 4, which is somewhat larger than expected |^. The global defect 



models normalized to the COBE results for the microwave background fluctuations 
probably yield somewhat too small perturbations in the dark matter. Nevertheless, 
the uncertainties concerning the bias factor and the nonlinear physics going into the 
calculation of the bias factor seem to me to leave room for doubts. It would be 
more convincing to rule out the the scenarios from the completely linear determina- 
tion of the microwave background fluctuations alone. So far, only the gravitational 
interaction of the radiation field with perturbations has been taken into account. 
To calculate 6T/T{6) on scales, 6 < 2°, which enter the horizon before recombina- 
tion, when baryons and photons still are a tightly coupled fluid, the baryon photon 
interaction has to be taken into account and the recombination process has to be 
modeled. For pure CDM without scalar field this calculation has been performed 
on different levels of accuracy [^, |39|, 0. For global defect induced fluctuations, 



intermediate and small scale anisotropies have only been approximated in the case 
when the universe is reionized at some early redshift, z > 100, and baryons and 
photons are coupled again via Thompson scattering. In this case, photon diffusion 
severely damps fluctuations on scales smaller than the horizon at z ~ 100, i.e., on 
all scales smaller than about 5° [^, [1^, ^. One of the missing pieces in the global 



defect scenarios is thus a detailed calculation of the microwave background spectrum 
on angular scales 6 < 2° or I > 100 for a non reionized universe. 

On the other hand, the CMB fluctuations, are not determined by the spectrum 
Q alone. The spectrum just yields the two point correlation function which deter- 
mines the fluctuations only if they are Gaussian distributed. In general the aim also 
yield non zero higher correlation functions. The skewness 5* and the kurtosis, K of 
the distribution of 2° x 2° pixels are found to be [^ 



S = -4 ±2.3 
K = 32 ± 29 . 

The deviation from Gaussian distribution is also shown in Fig. 8. These higher 
order correlations are an important mean to distinguish models with global defects 
from models with initial fluctuations from an inflationary epoch which usually yield 
Gaussian fluctuations. 



Figure 8. The pixel distribution of 1ST jT in a synthesized map of microwave back- 
ground fluctuations for the texture scenario of structure formation. The dashed hne 
shows a Gaussian with the same width and the same number of pixels. The negative 
skewness and the positive kurtosis are clearly visible, 
/.paragraph. 

5. CONCLUSIONS 

We have discussed gravitational interaction of global scalar flelds with matter and 
radiation. We have found that static global fleld conflgurations with vanishing po- 
tential energy do not affect slowly moving particles and do not redshift photons. 

Topological defects which form during phase transitions in the early universe 
can have important cosmological consequences. For gravitational interactions of the 
defects with the cosmic matter and radiation to be relevant, the defects must form 
due to a phase transition at GUT scale. In this case they may even seed the forma- 
tion of cosmological large scale structure. Even though it is not yet clear if defect 
induced structure formation scenarios do work out in detail, up today they remain 
an intriguing alternative to initial fluctuations from inflation since they also yield a 
scale invariant spectrum of perturbations. 



Acknowledgement: I'm grateful to my collaborator Zhi-Hong Zhou, 

who helped me in preparing some of the flgures for these proceedings. I have 
learned a lot from discussions with many of the participants, especially Robert Bran- 
denberger, Pedro Ferreira, Tom Kibble, Andrew Little, Leandros Perivolaropoulos, 
Tomislav Prokopec, Paul Shellard and Neil Turok. Finally, I want to express special 
thanks to Ann Davis, the organizer of this lively and stimulating meeting. 



REFERENCES 

[1] Langaker and Pi, Phys. Rev. Lett. 45, 1 (1980). 

[2] T.W.B. Kibble, Phys. Rep. 67, 183 (1980), see also these proceedings. 

[3] C. Nash and S. Sen, Topology and Geometry for Physicists, Academic Press, London 
(1983). 

[4] S. lyanaga and Y. Kawada, Encyclopedic Dictionary of Mathematics, Mathematical 
Society of Japan, MIT Press, Cambridge, MA (1977). 

[5] S.T. Hu, Homotopy Theory, Academic Press, London (1959). 

[6] J. Eells and L. Lemaire, Bull. London Math. Soc. 10, 1 (1978); 
ihid. 20, 385 (1980); 
R.T. Smith Proc. Amer. Math. Soc. 47, 229 (1975). 

[7] N. Turok and D.N. Spergel, Phys. Rev. Lett. 64, 2736 (1990). 

[8] G. Derrick, J. Math. Phys. 5, 1252 (1964). 

[9] L. Perivolaropoulos, Nucl. Phys. B375, 665 (1992). 

[10] S. Aminneborg, Nucl. Phys. B388, 521 (1992). 

[11] R.A. Leese and T. Prokopec, Phys. Rev. D44, 3749 (1991). 

[12] J. Borill, E. Copeland and A. Liddle, Phys. Lett. B258, 310 (1991). 

[13] T. Prokopec, A. Sornborger and R. Brandenberger, Phys. Rev. D45, 1971 (1992). 

[14] L. Perivolaropoulos, Phys. Rev. D46, 1858 (1992). 

[15] R. Durrer, Fund. Cosmic Physics 15, 209 (1994). 

[16] M. Barriola and A. Vilenkin, Phys. Rev. Lett. 63, 341 (1989). 

[17] R. Durrer, M. Heusler, P. Jetzer and N. Straumann, Nucl. Phys. B368, 527 (1992). 

[18] R. Durrer, Phys. Rev. D42, 2533 (1990). 

[19] M.J. Geller and J.P. Huchra, Science 246, 897 (1989). 

[20] G.F. Smoot et al., Astrophys. J. 396, LI (1992); 

E.L. Wright et al., Astrophys. J. 396, L13 (1992). 

[21] E.M. Lifshitz, Soviet J. Phys. (JETP) 10, 116 (1946). 

[22] E. Kolb and M. Turner, The Early Universe, Addison Wesley, Redwood City (1990). 

[23] R. Durrer and Z.-H. Zhou, Phys. Rev. Lett., submitted (1994). 

[24] T. Padmanabhan, Structure Formation in the Universe, Cambridge University Press 
(1993). 

[25] G. Efstathiou in: Physics of the Early Universe, Scottish Universities Summer 
Schools in Physics, ed. J. A. Peacock, A.F. Heavens and A.T. Davies (1990). 

[26] E.R. Harrison, Phys. Rev. Dl, 2726 (1970). 

[27] Ya. Zel'dovich, Mon. Not. Roy. ast. Soc. 160, Ip (1972). 

[28] P. Mezaros, Astron. Astrophys. 37, 225 (1974). 



[29] K.B. Fisher, M. Davis, M.A. Strauss, A. Yahil and J.P. Huchra, Astrophys. J. 402, 
42 (1993). 

[30] W. Press, D. Spergel and B. Ryden, Astrophys. J. 347, 590 (1989). 

[31] W.H. Press, B.P. Flannery S.A. Teucholsky and W.T. Vetterling, Numerical Recipes, 
Cambridge University Press (1990). 

[32] U. Pen, D.N. Spergel and N. Turok, Phys. Rev. D49, 692 (1994). 

[33] J. Borill, E.J. Copeland, A.R. Liddle, A. Stebbins and S. Veeraraghavan, Phys. Rev. 
D50, 2469 (1994), and references therein. 

[34] J. Borill, Phys. Rev. D50, 3676 (1994). 

[35] D. Bennett and S. Rhie Astrophys. J. 406, L7 (1993). 

[36] R. Durrer, A. Howard and Z.-H- Zhou, Phys. Rev. D49, 681 (1994). 

[37] N. Katz, L. Hernquist and D.H. Weinberg, Astrophys. J. 399, L109 (1992); 
R. Cen and J. Ostriker, ibid., L113 (1992). 

[38] N. Sugiyama and N. Gouda, Prog. Theor. Phys. 88, 803 (1992). 

[39] W. Hu and N. Sugiyama, Astrophys. J. Lett., submitted (1994). 

[40] U. Seljak, Astrophys. J. Lett., submitted (1994). 

[41] R. Durrer, J. Infrared Phys. Technol. 35, 83 (1994). 

[42] D. Coulson, P. Ferreira, D. Graham and N. Turok, Nature 368, 27 (1994). 



4- 



l_ / 

I 

I 

I 

I 
I 

I 

u 
I 



1 — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 
/\ 



\ 




J I L 



10 



1 - 



< 







-1 - 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




J I I I I I I I I I I I I I I I I 



-50 




T (t ~ 20) 



50 





1 


, ,V| 


iii| 


1 1 1 1 1 1 1 L 




:/ 






- 


0.1 


— 




\^ 


— 


CVJ 

CO 

CO 
Q. 0.01 


- 






II 1 1 1 1 III 


n nm 


1 


N 
1 1 1 


= 2 

M.I 


\ = 





10 



100 




0.001 



100 




0.001 



100 




0.001 



100 



im 




0.1 - 



0.001 



0.01 



k (Mpc-') 



P(k) 


- 1 1 


Mill 1 


1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 


II 1 1 1 1 


0.1 


- 




'^\ 




- 


0.01 


- 




^^ 


\ 


- 


D.OOl 


z— 






\ 


- 


DOOl 


z— 






^ 


\ 1 


10-5 


1 1 


1 


1 


1 






0.01 0.1 1 

k (Mpc-i) 



1000- 



X 

•1— I 




500- 



